Thermal Hall Conductivity of High Temperature Superconductors: 

Quantization and Scaling 
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Presented is the theory of thermal Hall conductivity K xy in the vortex state of high temperature 
cuprate superconductors in the clean limit. We show that limT->o Kxy/T is a staircase function of 
l/B with an envelope that scales as l/B. The relation to the experiments is discussed. 



Fifteen years after the discovery of high temperature 
superconductivity in cuprates the search for the "fi- 
nal" theory continues unabated. Nevertheless, certain 
important clues about the elements that such a theory 
must contain appear to be established. One of the most 
significant results of the past few years is that, in the su- 
perconducting state, cuprates seem to be well described 
by the familiar BCS-type formalism with Cooper pairs 
binding in the d-wave channel Q . This is a far from triv- 
ial observation in the materials known to exhibit strong 
correlations. Soon after the original discovery, Anderson 
H identified one of the key features of the new super- 
conductivity: the materials are basically two-dimensional 
with most of the action relevant for superconductivity 
taking place in the CuC>2 planes. The reduced dimen- 
sionality and the presence of nodal points on an otherwise 
gapped Fermi surface governs the low temperature prop- 
erties of the high-T c superconductors (HTS). Corrobo- 
rating evidence for the BCS character of the fermionic 
excitations comes from diverse spectroscopic and trans- 
port experiments ||, |6[ Q . 

Especially informative probe of the quasiparticle dy- 
namics is the thermal Hall conductivity, n xy g , particu- 
larly since phonons, which usually contribute to the lon- 
gitudinal thermal conductivity, k xx , do not contribute 
to K xy by the virtue of their being electrically neutral. 
In addition, from the theoretical point of view, in the 
leading "nodal" approximation, n xy vanishes || [lO| [ll[] 
and one must consider curvature terms. For these rea- 
sons a special opportunity emerges to scrutinize, both 
theoretically and experimentally, the extent to which the 
fermionic quasiparticles in HTS exhibit the BCS-type be- 
havior. 

Our starting point is a 2-dimensional Bogoliubov- 
de Gennes Hamiltonian in the vortex state 
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where A(x) is the vector potential associated with the 
uniform external magnetic field B, p is a momentum 
operator, and cr's are Pauli matrices operating in the 
Nambu space. The vortex phase fields satisfy V x 
V0(x) = 27rz^ i i5(x — Xj) with xj denoting the vortex 
positions and <5(x — Xj) being a 2D Dirac delta function. 
For a d x 2_ y 2-wave superconductor A(p) = ^r(p| — pV). 
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Following Franz and Tesanovic JO] (see also |Tl], |T^, 
H), we divide the vortices into two sets A and B, po- 
sitioned at {x^} and {xf } respectively (see Fig. |l|) 
and define two phase fields ^(x) and </> B (x) such that 
V x V</> Q (x) = 2ttz J2 t <K X - xf ), a = A,B. The Hamil- 
tonian (m is then transformed into |12 



H = U^H BdG U = <t 3 



(p - 03 v + a) 



2m* 



using the unitary operator 



U — g-iX+io-3* 



-cr 3 ^ + criA(p + a) 
(2) 

(3) 
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where \ = \{cj)A - 4>b) and $ = \{<j> A + <j) B ) = \<f>. 
In Eq. (|2|), v = \y<p — f A is the superfluid velocity 
and a = ^(V^ — ^4>b) is the Berry gauge field which 
imposes the condition that a quasiparticle wavepacket 
encircling an |p vortex must pick up an overall minus 
sign. 

In addition to keeping the wavefunctions single- valued, 
the advantage of this "bipartite" singular gauge trans- 
formation is that for a periodic arrangement of vor- 
tices (Abrikosov vortex lattice) the resulting Hamiltonian 
(2) is itself periodic. This is not true of Hamiltonian 
(1). Therefore, the eigenstates of (||) can be expressed 
through the familiar Bloch waves ?Ank(x) — e lk x u„k(x) 
labeled by one discrete quantum number, the band in- 
dex n, and two continuous quantum numbers, the vortex 
crystal momenta k x , k y that range in the first Brillouin 
zone. Note that u„k(x) is not periodic with the period- 
icity of the vortex lattice, but with the periodicity of the 
vortex sublattice containing a unit of electronic flux — 
(see Fig. 0). 

By extracting the plane wave part of the Bloch states 
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FIG. 1: Example of A and B sublattices for the square vortex 
arrangement. 



we can define a Hamiltonian 



W(k) = a 3 



(p + k cr 3 v + &f 
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- <7 3 fl + CTiA(p + k 



a) 

(4) 

where H(k) = e lkx 7ie lkx . We now wish to point out 
an important feature of 7i(k) in the case when the vortex 
lattice has an inversion symmetry. As shown in Ref. fl4f| 
an operator V = ia^CI, where C represents charge conju- 
gation and X takes x into — x, with the following property 



V- l H(K)V = -H(k); (u nk \V\u nls ) = 0. 



(5) 



This means that for every energy band e„k there is an- 
other e m k such that e„k = — e m k, i-e. the band spectrum 
has reflection symmetry about the Fermi level and if a 
band were to cross the Fermi level (e nq = for some 
q) then there would have to be a twofold degeneracy at 
q. However, by the von Neumann- Wigner "non-crossing" 
theorem [|l5j , the degeneracy unrelated to symmetry can- 
not happen for a time reversal breaking Hamiltonian 
which depends only on two continuous parameters. Thus, 
the quasiparticle spectrum of the d x 2_ y 2 superconductor 
in the vortex state is gapped by a (direct) gap A m <C Ao 
unless a third parameter is fine tuned, e.g. the band 
crossing at the Fermi level can be achieved by fine-tuning 
the chemical potential /i. 

If we rescale all coordinates by the magnetic length 

Zx, 



I = y/c/eB (see Fig. 1), x 
W(k,Z,/x,A ) - 



then 



(6) 



where / is a universal (operator) function. Consequently, 
the small band gap A m at the Fermi level, which is equal 
to the minimum of the lowest positive eigen-energy of 
7i(k, Z,/i) with respect to k, is also a universal function 
that can be written as 



A m (B,^A )=Bx F 



fj, 



Ao 

PF 
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FIG. 2: Scaling of the field induced quasiparticle gap A m 
(open circles) with the magnetic field and the chemical poten- 
tial (measured from the bottom of the tight-binding band). 
The quantized values of the <7^ m are also shown (solid line). 
Magnetic length, I, is in units of the tight-binding lattice con- 
stant a and energies are in units of the hoping energy t; the 
results are for a square vortex lattice with Ao = O.lt. The 
dashed line is the envelope scaling function. 



Thus, in order to determine the dependence of A m on the 
magnetic field B, all we need to know is the dependence of 
A m on [i for some fixed B and , and then use the scal- 
ing relation (Q). This dependence is determined numeri- 
cally by diagonalizing a tight-binding version [ fl4| of the 
Hamiltonian (|j) close to the bottom of the tight-binding 
band where we expect to recover the continuum theory 
(JI|). The numerical calculations confirm the scaling (0) 
which gives us confidence in relating the continuum re- 
sults to the tight-binding approximation. As shown in 
Fig. ||, for a fixed value of^=a£i = 10, where t is the 
tight-binding hopping constant, -gA m (-^) is a saw-tooth 
function with a period of approximately An. 

In Ref. jfrj it was demonstrated that the gapped quasi- 
particle spectrum leads to the quantization of the spin 
Hall conductivity cr^ m (see also Ref. [[ll]), which is in 
turn related to the thermal Hall conductivity, n xy , via 
the Wiedemann-Franz law 
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Because the quantization of <T x v y m is topological in ori- 
gin fll4| |, any continuous change in a parameter of the 
Hamiltonian results in the same quantized value of <J x v y m 
unless the quasiparticle gap collapses. Therefore, for a 
fixed Dirac anisotropy an, by virtue of Eq. @, the de- 
pendence of cr s x p y n on B and n can be expressed entirely 
as a dependence of a x v y m on jjA m (B, fi), i.e. 



oZ n {B,n)=G 



(9) 



where G is a universal function of its argument that in 
addition depends on the geometry of the vortex lattice. 
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FIG. 3: Thermal Hall conductivity data of Zhang et. al. 
for B > 4T plotted vs. B" 1 . 



For a fixed (large) value of B the numerical computa- 
tion of <J x p y n can be readily accomplished. The results 
for a square vortex lattice and ao — 10 are displayed in 
Fig. H where it is seen that o x v y n is a stair-case func- 
tion of 4| which starts at zero at empty filling and then 
jumps by a multiple of 2-$=. The band crossings, which 
are responsible for the changes in a s x p y n , occur at symme- 
try points in the Brillouin zone, and by the symmetry of 
the square vortex lattice, the steps are bound to come in 
even integers Q . 

Thus we arrive at the main result of this Letter shown 



in Fig. 
by Eq. 



the envelope of the o- 8 x p y n (B) scales as B 



and 



with non- vanishing longitudinal quasi- 
Even in the presence of 



i.e. k xv - p* x+piv 
particle thermal resistivity p x 
weak disorder, we still expect p xy oc B and p xx to be ap- 
proximately field independent. Therefore, the transverse 
thermal conductivity in the presence of disorder should 
vanish in the limit of vanishing B. We must admit that 
the origin of the residual thermal conductivity at = 0, 
evident in Fig. ||, remains mysterious to us and we can 
only speculate that it is not of quasiparticle origin. 

In summary, based on a general analysis of the full 
non-linearized Bogoliubov-de Gennes equations, we ar- 
gued that magnetic field B will induce a small gap A m 
in the quasiparticle spectrum which is a non-monotonic 
function of the magnetic field. Rather, A m vanishes at 
some special values of fi/B which depend on the details 
of the vortex lattice and the Dirac cone anisotropy (Fig. 
H). The dependence of n xy on B is entirely due to this 
non-monotonic behavior of A m , since n xy /T is a topolog- 
ical quantity and therefore can change only at the band 
crossings ]T4| . The topological quantization of n xy was 
also discussed in the Ref. |ll[] . However, we wish to point 
out a significant qualitative difference between the results 
presented here and those of Ref. jllj]: as we argued be- 
fore quantized staircase function of 1 / B with an 
envelope that itself scales as 1/B, while Ref. [[n] states 
can take on only one of the three allowed 
values 0, ±2 and is otherwise field- independent. Conse- 
quently, in a typical experimental situation, the explicit 
values for n xy obtained from these two works will differ 
significantly. Still, there is an agreement on the fact that 
since the quasiparticle spectrum is gapped n xy does not 
obey simple Simon-Lee scaling, at least when the tem- 
perature is comparable to the quasiparticle gap. 

We thank Marcel Franz, Ashot Melikyan and Ady 
Stern for comments and discussions. This work was sup- 
ported in part by NSF grant DMR00-94981. 
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This expression is reminiscent of the electrical Hall con- 
ductivity in the 2D electron gas. 

Recently, measurements of n xy were conducted by 
Ong's group M on YBCO samples with a very long mean 
free path. These experiments were carried out over a 
wide range of magnetic fields (up to ~ 14T) and at tem- 
peratures from T ~ 12. 5K to above the superconducting 
transition T c ~ 90-fT. Unfortunately, the experiments 
are resolution limited below 12.5K as signal becomes too 
weak. In Fig. || we re-analyze the data of Ref. Q . It is 
seen that for the magnetic fields above 4T, n xy oc C+B^ 1 
which is, apart from the finite intercept, in agreement 
with our theory. At fields B < 4T, n xy decreases with 
decreasing B and vanishes at B = jsj. At low mag- 
netic fields -i divergence of K xy is cut-off by disorder 
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